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Abstract. Motivated by field-theoretic predictions we investigate the stable excitations that exist in two 
characteristic gapped phases of a spin-1 model with Ising-like and single-ion anisotropies. The sine-Gordon 
theory indicates a region close to the phase boundary where a stable breather exists besides the stable 
particles, that form the Haldane triplet at the Heisenberg isotropic point. The numerical data, obtained by 
means of the Density Matrix Renormalization Group, confirm this picture in the so-called large-D phase 
for which we give also a quantitative analysis of the bound states using standard perturbation theory. 
However, the situation turns out to be considerably more intricate in the Haldane phase where, to the best 
of our data, we do not observe stable breathers contrarily to what could be expected from the sine-Gordon 
model, but rather only the three modes predicted by a novel anisotropic extension of the Non-Linear Sigma 
Model studied here by means of a saddle-point approximation. 

PACS. 75.10.Pq Spin chain models - 11. 10. St Bound and unstable states; Bethe-Salpeter equations - 
ll.10.Kk Field theories in dimensions other than four 



1 Introduction 

One-dimensional spin systems have been largely studied 
since Haldane proposed his conjecture. He suggested 
that while half-odd integer spin chains should always be 
gapless, integer-spin ones should manifest a gap. The con- 
jecture is supported by a field-theoretic analysis [2] of 
the lattice Hamiltonians, a typical approach for low-di- 
mensional quantum systems that admit a continuum-limit 
counterpart in (D + 1) dimensions. The low-energy spec- 
trum of the spin Hamiltonian can be interpreted in the lan- 
guage of (quasi)particles: a finite energy gap corresponds 
to a massive particle at rest, and the dispersion relation 
of the excitation for small momenta can be read as a rel- 
ativistic (on shell) energy-momentum relation. The field- 
theoretic approach proved to be extremely powerful to ex- 
plain experimental data. For instance, recent observations 
on spin-1/2 compounds with Dzyaloshinskii-Moriya inter- 
action in external field [SHE], confirm the appearence of 
effective particles, namely solitons and breathers, as pre- 
dicted by the Sine-Gordon Model (SGM) that describes 
the low energy spectrum of the chain. 

In this paper we discuss the possibility of observing 
solitons and possibly breathers in spin-1 chains with in- 
ternal anisotropies. In particular, it is important to under- 



stand the role of the interactions between the particles of 
the continuum theory to see how their bound and scatter- 
ing states manifest in the energy spectrum of the lattice 
model. The effective interactions will depend on the values 
of parameters of the spin Hamiltonian and it is possible 
that certain two- (or more) particle states that arc stable 
in a region of parameter space loose their stability when 
the parameters are continuously changed. In other terms 
the formerly stable state acquires a rest mass equal or 
larger than the sum of the rest masses of its constituents 
and, if there are no special selection rules, the particle 
decays into the continuum part of the spectrum. 

In this paper we study such a scenario for the following 
spin-1 anisotropic model: 



H = 



qy qy 



\S*Sz +1 +D(S*) 2 }, 
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where A and D parametrize the Ising-like and single-ion 
anisotropies respectively. This Hamiltonian is known to 
be described, in the continuum low-energy limit, by a 
Non-Linear Sigma Model (NLcM) in the vicinity of the 
isotropic antifcrromagnctic Heisenberg point (A = 1, D = 
0) HIE]- and the SGM in the neighborhood of the critical 
line separating two gapped phases, namely the large-D 
and the Haldane ones 0- These phases will be defined 
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in section [21 where we will also recall the mapping onto 
the SGM and present a novel extension of the NLctM 
for anisotropic integer-spin models which encompasses the 
usual isotropic NLerM that lies at the basis of the Haldane 
conjecture. In general we find a triplet of excitations and 
no other bound states. On the contrary, from a quantita- 
tive analysis of the SGM we expect that in some region 
of the Haldane and large- D phases at least an additional 
bound state (a breather) should appear. We proceed in 
section to a numerical check of these expectations using 
the multi-target Density Matrix Renormalization Group 
(DMRG) technique that allows us to handle chains of up 
to 100 sites and extract various excited states of the spec- 
trum. The data in the large-D phase confirm the existence 
of a breather in the region predicted by the SGM. On the 
other hand, we do not find any stable breather in the 
Haldane phase. In section 2] we will sketch also a simple 
perturbative argument that provides a quantitative inter- 
pretation of the data. Finally, in section[3]we will comment 
on our results and draw some general conclusions. 




A 



2 Quantum field theories for the A — D model 

Among integer spin models, the anisotropic S = 1 chain 
with its rich phase diagram occupies a relevant position. 
At the isotropic point, (A = 1, D = 0) the 0(3) Hciscn- 
berg model is recovered. The full phase diagram consists 
of six different phases (see [7| for a recent numerical de- 
termination). We will focus our attention on the A > 
half-plane where, apart from the transition lines, all the 
phases show a nonzero energy gap above the ground state 
(GS). In fig.rj]we report the phase diagram in this range of 
parameters. Usually, two different types of order parame- 
ters are used to characterize these phases: the Neel order 
parameters (NOP): 



0%= lim (-1)*-'(S?S;); a = x,y,z (2) 
and the string order parameters (SOP): 



lim (Sfe 1 

-jj-»oo 



: ^=^ S *S?); a = x,y,z, (3) 



first introduced by den Nijs and Rommelse For D>1 
and D > A the system is in the so-called large-D phase, 
with a unique GS that does not break the above symme- 
try: Og = 0% = for all a. Decreasing D, for A < 3, we 
have a transition into the Haldane phase. This phase is 
characterized by the non-vanishing of all the components 
of the SOP's Og 0, meaning that the Z2 x Z2 symmetry 
is fully broken, and by 0^ = Va. Inside the Haldane 
phase, the gap of the first excited state, though always 
different from zero, belongs to two different spin sectors 
depending on the parameters A and D. It belongs to the 
5f ot = ±1 sector for values of A and D above the so-called 
degeneracy line (see fig.^l. This line, which passes through 
the 0(3) isotropic point, divides the Haldane phase into 



Fig. 1. Phase diagram of the model Q obtained numerically 
via DMRG. The two lines passing through the isotropic point 
are the loci of points where the transverse and longitudinal 
gaps coincide. Circles are DMRG data whereas squares are the 
results of the theoretical approach. The vertical lines across 
the Haldane-large-D transitions join the points studied numer- 
ically (see section |3- 



two sub-phases (2). Below this line the first excited state 
has S* ot = 0. The GS is always unique and belongs to the 



Of , 



sector. On the right of the Haldane phase we 



find the twofold degenerate Neel phase. Here the Z2 x Z2 
symmetry is only partly broken as as proved by the fact 
that 0% = = for a = x,y but O z N , O z s ^ 0. The 
Haldanc-largc-D (H-D) and Haldanc-Ncel (H-N) transi- 
tion lines mark second order transitions; the former is de- 
scribed by a c = 1 Conformal Field Theory (CFT) while 
the latter by a c = 1/2 CFT. They merge at the tricritical 
point (A = 2.90, D = 3.20) [5] where the Haldane phase 
disappears and the large- -D-Neel transition becomes first 
order. 

Near the isotropic point the Hamiltonian (TJJ can be 
mapped onto an anisotropic version of the NLerM using 
Haldane's ansatz FJ] to decompose the spin coherent states 



field Qj{t) = nj(T)(-iyyi- 



*3\' J — "3V 1 J\ V S 2 

form, £j (r) , and a staggered part rij (r) 
satisfying the constraint |njj 2 + r 



i(r) 



into a uni- 
[n 3 t(T),n 0Z ( T )], 
1. After integration 
of the uniform field, one obtains the following effective 
Lagrangian density: 



C 



Vz 

2ff z 



{d x n z ) +^{d T n z ) + ^, z ni 



+ A„t, (4) 



where the interaction term has a simple expression only 
in the regions \n z \ 1 and |njj <C 1. In the first region 
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the interaction term has the form 



C mt =nl d_i\d T nj_\ + d z (d T n z ) 2 , 



In a previous work jHj it was argued that near the 
Gaussian c = 1 line the perturbed 0(2) model can written 
in terms of a sine-Gordon model: 



whereas, in the region |nj 
the form 



<C 1, the interaction term has 



&nt = Inj 



d±\d T n_ L \ +d z {d T n z y 



Here the constants g±( z ), ^±(2)7 ^_l(z)j A*z are different func- 
tions of the microscopic parameters A and D whose ex- 
plicit expressions depend on the region one considers. 

As long as one deals with integer spins the topological 
term arising from the combination of the Berry phases on 
each site reduces to a multiple of 27r and therefore can 
be safely omitted. More details on the mapping into this 
anisotropic variant of the NLctM can be found in ref. ^U] ■ 
The constraint of unit norm makes the theory hard to 
treat. In both regions, \n z \ <C 1 and |nj_| € 1, we limit 
ourselves to a mean-field solution that can be worked out 
using the saddle-point approximation on eq. The con- 
straint is taken into account by introducing a uniform La- 
grange multiplier r\. We refer the reader to ref. ^Uj for de- 
tails and only mention that one arrives at a couple of self- 
consistent equations for the longitudinal and transverse 
gaps AE Z and AE±, that play the role of the masses of 
the particles of an "anisotropic Haldanc triplet" . Of course 
AE Z = AE± at the isotropic Heisenberg point. However, 
when A and D are varied the transverse and longitudi- 
nal channels split. Nonetheless, guided by the numerical 
work of [§], we can search for a specific line in parameter 
space where the two gaps remain degenerate (in the ther- 
modynamic limit) even if the lattice theory has no 0(3) 
symmetry. Interestingly, the line found through the self- 
consistent equations is in good agreement with the one 
found with the DMRG as depicted in fig. ^ 

Let us now see a possible connection with the SGM. 
Neglecting the longitudinal field (i.e. setting n z = 0) and 



writing 



,M 



then an 0(2) NLctM is recovered: 



C 



(d x ey 



4 (drOf 



(5) 



with = 9/g, g = ^/(2 + 2\ + D) and v = g (for 5 = 1) 
. This is exactly a free Gaussian model with a bosonic 
field compactified along a circle of radius lf-^/g. We know 
that this model describes a CFT with central charge c = 1 
and with primary fields having scaling dimensions: 



AK 



n 2 K 



(6) 



with K = ir/g and m, n £ Z. 



1 In ref. |S] in all the expressions after eq. (6) D has to be 
replaced with D/2. This replacement affects the theoretical val- 
ues of K and of the scaling dimensions, but not their numerical 
estimates. 



C 



v (d x $f + - (d T $f 

V 



Vfl 



cos (V4nK$) (7) 



where <P is a compactified field dual to 0, a is the lattice 
spacing and /2 is a coupling constant that should vanish 
exactly along the H-D transition line. K is the basic pa- 
rameter that allows to compute all the scaling dimensions, 
and consequently all the critical exponents, that along 
the c = 1 line acquire nonuniversal values depending on 
the actual values of A and D: it ranges from the value of 
K = 2 at the Berezinskii-Kosterlitz-Thouless (BKT) point 
(A = 0, D ~ 0.4) to the value K = 1/2 at the tricritical 
point, becoming K — 1 at the so-called free Dirac point, 
A ~ 2. The SGM has a spectrum consisting of a soliton, an 
anti-soliton and their bound states, breathers, the number 
of which is different from zero only for K < 1 |llj . This 
means that, besides the excitations of the c = 1/2 theory 
(belonging to S* ot = 0) , in the region of the phase diagram 
that remains on the left side of the free Dirac point (i.e. ap- 
proximately for < A < 2, for which 1 < K < 2) the par- 
ticle content of the theory should consist only of a soliton 
and an anti-soliton. When K < 1 (i.e. 2 < A < 2.90), also 
a breather should appear (and possibly a second breather 
for K > 2/3). If this is proved to be correct, then it would 
mean that the Haldane and the large- D phases could be 
eventually classified into subphases according to the num- 
ber of stable particles. In order to check this hypothesis, 
we have numerically analyzed the spectrum of the model 
JIJ on the A = 1 and A = 2.59 lines for values of D crossing 
the Haldane-large-D transition curve. 



3 Numerical investigation 

3.1 K > 1: No breathers 

For A = 1 the critical point on the c = 1 line has been 
previously located at D c = 0.99 with a parameter K = 
1.328 ± 0.004 UJ; then we are in the sector of the Haldane 
phase that should be characterized only by a soliton and 
an anti-soliton. From CFT we know that a free Gaussian 
theory for a field compactified along a circle with radius 
K/tt has primary fields of scaling dimensions given by 
equation JHJ where m plays the role of S' ot . The energies 
of the excited states are related to the scaling dimensions 
by the formula: 



AE mn E mn Ei 



2ttv 

: = — (d mn + r + f) , r, f € N. (8) 



The spectrum of the scaling dimension is reported in ta- 
ble 2 of ref. 0. In the language of the SGM, labelling 
the states as (m,n,r,r), the first excited states are given 
by (±1,0,0,0) and belong to the S* t z ot = ±1 spin sectors. 
They correspond to the soliton and the anti-soliton, while 
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0.02 0.03 

1/L 



0.05 



Fig. 2. Plots of the low-lying energy gaps versus 1/L for D = 
0.90 A = 1. The states are labeled with the quantum numbers 
of the CFT (see text). 



the first breather comes from the doublet (0, ±1, 0, 0) that 
splits into two different states as soon as one moves away 
from criticality. Operatively, a breather is a singlet in the 



spin sector. To be a stable state, the ratio between 



its rest energy and that of the soliton (measured respect 
to the ground state) R = AE^j AE S must be smaller than 
two. 

The numerical analysis has been performed with a 
DMRG algorithm that exploits the so-called thick-restart 
Lanczos method (see ^2] for a recent review and ^2] for 
our implementation). We adopt periodic boundary condi- 
tions (PBC) to avoid complicancies arising from edge ef- 
fects (midgap states and surface contributions see e.g. |14|). 
and perform up to five finite-system iterations in order 
to reduce the uncertainty on the energies to the order of 
magnitude of the truncation error ^H]- Using M — 400 
DMRG states the latter is O(10 -4 ) or better, the worst 
cases being the ones with f target states. 

On the A = f line we have studied points close the 
Haldanc-largc--D transition as well as points inside the 
gapped Haldane and large- D phases. For D = 0.85, 0.90, 
0.95, L06, 1.11, we are still in a quasi-critical region for 
which the correlation length is still larger than the size of 
the system. This is confirmed by the fact that the gaps 
scale substantially linearly in 1/L as predicted by eq. © , 
and we can label the states with the quantum numbers of 
the CFT at the critical point. A typical situation is shown 
in fig. |2 and one sees that the states corresponding to 
(0, ±1,0,0) are still degenerate. Following the evolution 
of this doublet for values of D farther from the critical 
line, we can identify the first breather as the lowest of the 
two splitted states. 

Moving downward with D well into the Haldane phase 
up to D = 0.3, 0.6, in the S% ot = sector we observe be- 
sides a doublet state that we identify with the secondary 




0.2 



(+1,0,0,0) 
(+2,0,0,0) 
(0,+l,0,0) 
© — © (0,-1,0,0) 
b-d (0,0,1/0,0/1) 
x^< c = 1/2 
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Fig. 3. Plots of the low-lying energy gaps versus L for D — 
0.3 at A = 1. States are still labeled as in fig. |5] except that 
here the states (0, ±1, 0, 0) have split and grey circles represent 
the highest of the two. It is already possible to see the first 
excitation of the c = 1/2 theory (crosses); this state appears 
incomplete because for L > 50 it lies outside the set of the 
DMRG-targeted states. 



states (0, 0, 1/0, 0/1), two clearly separated states that are 
the evolution of the (0, ±1,0,0) doublet. This feature is 
evident in fig. reporting the case D = 0.3. At these 
values of D the correlation length is sufficiently small, 
and we extrapolated the asymptotic values of the energy 
gap for both the soliton and the breather using the semi- 
phenomenological formula |16| 



AE{L)=AE(po) + j- expf-^j 



(9) 



where the fitting parameter £' is expected to be propor- 
tional to the actual correlation length. We verified that 
eq. © , fits better the data than other functions with dif- 
ferent exponents of 1/L. The extrapolated values for the 
breather gaps lay slightly above the continuum threshold, 
specifically we obtain R = 2.048 at D = 0.6 and R = 2.043 
at D = 0.3. In addition we verified that also the lowest 
state in the S* ot = 2 sector remains in the continuum in 
the thermodynamic limit. 

We pushed the analysis a little further for D > 1.5. As 
the correlation length is now sufficiently small we trust the 
data at L = 100 without the help of any extrapolations. 
The result is that also in this regime R > 2 so that there 
are no stable breathers. In agreement with what observed 
in |17U18j . where the finite-size spectrum of the SGM is 
studied, the first excitation in the sector 5* ot = lies 
always above the soliton. 

In summary, as far as the K > 1 region is concerned, 
the numerical investigation allows to conclude that there 
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d u * 1 [ x degeneracy] 


(m, n, r, r) 


[xl] 


(0,0,0,0) 


0.294 ±0.003 [x2] 


(±1,0,0,0) 


0.85 ±0.01 [x2] 


(0,±1,0,0) 


1 [x2] 


(0,0,1/0,0/1) 


1.18 ±0.01 [x2] 


(±2,0,0,0) 


1.29 ±0.02 [x4] 


(±1,0,1/0,0/1) 



Table 1. Spectrum of scaling dimensions at the point (A = 
2.59,D C = 2.30), obtained from formula © where K = 0.85 ± 
0.01 El. 



are no soliton bound states both in the Haldane and in 
the large- D phases. 



3.2 K < 1: Emergence of the breather 

Let's analyze now what happens on the A = 2.59 line, 
for which we have previously checked that D c = 2.30 and 
K = 0.85 (Bj. If we believe in the SGM as a faithful con- 
tinuum theory of the low energy part of the model JQl 
we expect the spectrum to present a stable breather both 
in the Haldane and in the large-D phases. According to 
the estimated K the order of the CFT energy levels is the 
one reported in tabled where we have shown the scaling 
dimensions 

d CFT = 

As in the case K > 1, we are interested in the soli- 
ton gap originating from the (±1, 0, 0, 0) states and in the 
breather coming from the (0, ±1,0,0) states. 

The correlation length now decreases rather rapidly 
with \D — D c \ so that we see the off-critical regime for D 
quite close to the critical point for the system sizes at our 
disposal. 

Inside the large- D phase we selected eight points with 
D ranging from 2.39 to 3. We have checked that the corre- 
lation length is indeed so small that we can use the data at 
L = 100 without any finite-size scaling. The identification 
of the soliton and the breather states out of the DMRG 
spectrum is rather direct. Apart from a single point quite 
close to the critical line, the gap ratio R is always less than 
two, as shown in fig. 0] This confirms the existence of an 
additional stable particle that we identify as the breather. 
Considering the non-monotonic behavior of the function 
R(D), we can speculate whether the breather looses its 
stability for very large values of D. In the next section we 
will study this point by means of an analytical approach. 

In the Haldane phase, which is very narrow here, we 
studied the points D = 2.2222, 2.235, 2.25. The identifica- 
tion of the soliton and the breather is now complicated by 
the appearance of states which do not originate from the 
c = 1 CFT. A typical numerical spectrum up to L = 100 is 
shown in fig. El Following the states in the S* ot = sector, 
labelled by crosses, for smaller values of D we observe that 
their gaps decrease and eventually vanish at the c = 1/2 
transition line whereas all the other states remain massive. 
For this reason we identify these states as excitations com- 
ing from the c = 1/2 CFT. The doublet (triangles down) 




D 



Fig. 4. Plot of the ratio R between the breather and the soliton 
energy gaps at L — 100 as a function of D for A = 2.59. 
The dashed horizontal line represents the continuum threshold 
while the dotted vertical line marks the critical point. 




Q| 1 1 1 1 1 1 1 

20 40 60 80 100 

L 

Fig. 5. Typical low-lying spectrum in the Haldane phase for 
K < 1 (A = 2.59, D = 2.235). 



cannot be identified with the breather being twofold de- 
generate; it is probably a descendent (r = 0/1, f = 1/0) 
of the first c = 1/2 state. What we find is that, both 
at D = 2.25 and at D = 2.235, the breather lies above 
the two-soliton gap continuum, i.e. R > 2. The point 
D = 2.2222 lies on the "degeneracy-line" where the soli- 
ton gap is degenerate with the level originating from the 
first excitation of the c = 1/2 theory. Now many states 
fall below the breather so that the latter can no longer be 
seen within the numerically available levels. 
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State 


Energy 


Degeneracy 


Name 


|0i---0i) 





1 


GS 


|0i ••■+,■• ■ L > 
1 < 3 < L 


D 


L 


Soliton 


|0 X hi k---0 L ) 

l<j^k<L 


2D 


L(L-l) 


Breather 



Table 2. Low energy states of the unperturbed Hamiltonian 



of the interaction the soliton states may form bound state 
if their energy lies below the two-particle continuum. At 
first perturbative order the breather energy can be cal- 
culated using the Bethe-Ansatz approach; details can be 
found in Appendix A. The resulting breather states arc 
indexed by the center-of-mass momentum Q. Acceptable 
states are those for which cos (Q/2) < A/2. This condition 
defines a characteristic wavenumber Q* above which the 
breather form a bound state with energy 



4 Perturbation theory at large- D 

From the data reported in fig. 0] one is led to separate 
the large- D phase into two regions depending on the sta- 
bility of the breather. However, the location of the curve 
A s tab(-D) separating the two regions is still vague: it starts 
from the c = 1 H-D line at the point where K — 1, close to 
A ~ 2, but we have no precise indications for larger values 
of D. In this region a perturbative analysis can be carried 
on. From the results quoted in |19j it is known that bound 
states (at zero total momentum) do not exist for A = 1. 
Here we show that the presence of an Ising-like anisotropy 
A allows for a stabilization of breather states. 

We start by treating the D-term as unperturbed Hamil- 
tonian 



= 2D 



H = Dj2(Si 

3=1 

and the rest as small perturbation 

L 



(10) 



.7=1 L ' 



A - - cos (Q/2) 
A 



O (D- 1 ) , for |Q| > Q* 



(11) 

Note that for Q = it the bound state is always below the 

(2) 

continuum with an energy shift SE^ = —A. However with 
the DMRG we inspect the low-lying part of the spectrum, 
and in order to find a stable bound state, that will be 
identified with the breather, we must examine the mini- 
mum of the continuum at Q = 0. In this case the bound 
states emerges from the continuum only for A > 2 with an 

(2) 

energy shift 5Eq = —A — 4/ A. The mechanism is illus- 
trated in fig. El for some values of A. Interestingly enough, 
the mass ratio R = AE\,(oo) / AE s (po) now reads: 



R 



E, 



(2) 



i (A -2y 

A (D — 2) 



(12) 



This equation is valid for A > 2, meaning that the sta- 
bility boundary is A s tat>(-D) = 2 independent of D. This 
result is consistent with the stability criterion provided 
by the SGM: K = 1 just for A - 2 0. In addition the 
ratio converges to R = 2 from below in agreement with 
fig- El The first-order approximation is also quantitatively 
good in comparison to the DMRG values as shown by the 
example in fig. 



Both Hq and V commute separately with the total z- 
component of the spin, S* ot , and with the total "spin- flip" 
(or "time- inversion" ) operator T = cxp(i7r5f ot ). However, 
[5 t 2 ot , T] ^ and one can specify both quantum numbers 
for the energy eigenstates only in the sector S£ ot = 0. 
The GS falls precisely in this sector and here we expect 
to find the breather. The soliton and the anti-soliton have 
S* ot = ±1, but we can restrict to the soliton case 5f ot = 1 
since the anti-soliton state is simply obtained by applying 
T. 

In table |2 we show the GS together with the first ex- 
cited states for the unperturbed Hamiltonian 1)1 ()() . In this 
picture the "solitons" are completely degenerate and non- 
interacting so that the (soliton-antisoliton) "breather's" 
energy is exactly twice that of the 1-particle states. As 
we switch on the interaction V, the GS energy becomes 
E^ = —L/D + O (D~ 2 ) up to second order in perturba- 
tion theory. The first excitations are no longer degenerate; 
they rather form a band, labelled by the lattice momentum 
q, E ( q ] = D + 2 cos (q) + O (Z)" 1 ) . Up to first order per- 
turbation theory these 1-particle states form a continuum 
of two-particle scattering states whose upper and lower 
edges arc given by the dashed lines in fig. As a result 



5 Discussion of the results and conclusions 

Motivated by field-theoretic predictions, in this paper we 
have investigated the stability of massive excitations (soli- 
ton, ant-isoliton and breathers) that are present in the 
large- D and Haldane phases of the anisotropic spin-1 Heisen- 
berg chain with both Ising-likc and single-ion anisotropy. 

In excellent agreement with the SGM that is supposed 
to describe the continum limit of the lattice model in the 
vicinity of the H-D critical line, we have provided both 
numerical and analytical results supporting the existence 
of a stable breather in the large-D phase of the model (QJ 
in the region with A > 2. 

On the contrary, as far as the Haldane phase is con- 
cerned, our best estimates do not show other stable states 
in addition to the three particles that form the Haldane 
triplet at the isotropic Heisenberg point. This is at vari- 
ance with what we would expect from the SGM, indicating 
that the latter is probably not sufficient to describe the ef- 
fective interaction when there is an underlying topological 
order as the one found in the Haldane phase (20j. Indeed 
in this phase we may expect that we cannot neglect the 
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Fig. 6. Bound states (full lines) and continuum (region en- 
closed by dashed lines) for the values of A indicated on the 
right. Full dots mark the wavenumbers Q* where the bound 
state starts to detach from the continuum when A < 2. 
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Fig. 7. Results of first-order perturbation theory vs DMRG 
data in the large-D phase at A = 2.59. 



effects of an interaction between the states coming from 
the c — 1 CFT describing the H-D critical line and those 
pertaining to the c = 1/2 CFT describing the H-I transi- 
tion. 

It is interesting to speculate whether this excitation 
can be observed in real materials. The so called NENC 
has the advantage of having a large single-ion anisotropy 
D ~ 7.5 but A = 1 in this compound. On the other 
hand large Ising anisotropics are present in some spin- 
1/2 materials (e.g. CSC0CI3 and CsCoBr 3 ; see an d 
references therein). Interestingly enough, differently from 



existing experiments the mechanism proposed here to ob- 
serve breather states does not require the presence of ex- 
ternal magnetic fields. 
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Appendix A 

In this appendix we show how to evaluate the first order 
correction of the second energy level in S* ot = 0. The 
unperturbed states are 



45) 



|0i 



/.•••"/. • J^k 



(2) (2) 

and we need the eigenvalues of the matrix (u\ i„JV]u!- Q 
with j ^ k and I ^= m. We write the -unnormalized- 
eigenfunction as 



j,k 



The states can be classified according to the eigenvalue of 
T: 



±1 



fkj — T fjk- 



In terms of the coefficients fjj. the eigenvalue equation 
reads: 

(1 - Sjk) [fj+l,k + fj,k-X + fj-lM + fj,k+l] 

- \f jk (5 jtk+1 + S jtk -i) = SE^frk. (13) 

Note that we have put the term (1 — Sjk) in the left- 
hand side for compatibility with the requirement fjj = 0. 
Now, eq. (|13(l is an eigenvalue equation in a L(L — 1)- 
dimensional space. It is actually much more simple to ex- 
press the amplitudes using relative and center-of-mass co- 
ordinates: r = j — k and p = (j + k)/2. From now on we 
will follow closely the approach used in Mattis' book |22I 
(sec. 5.3) for the bound states of magnons in Hciscnberg 
ferromagnets. Consider the wave function: 

f ik = cxp (iQp) F(r), F(r) = ^ f(p) exp (ipr) , (14) 



where Q and p denote the total and relative momenta 
respectively, both ranging in (— n, n]. The eigenvalue of 
T reflects directly in the parity of the Fourier transform: 
f{~P) = t /(p)i while the constraint fjj = becomes 
= J2 P f(p) = °- Whcn d is Pegged into 
a common factor exp (iQp) can be simplified from both 
members and we are left with the simpler equation: 



(1 — S r fi) cos 



Q 



[F(r + 1) + F(r - 1)] 



XF(r)(S ril +S 1 ,. 1 )=SE^F(r), (15) 
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that appears as a finite-differences Schrodinger equation 
with a (5-potential at r = ±1. From now on we fix r = 1; 
the procedure for r = — 1 is analogous and the conclu- 
sions are exactely the same. In momentum space eq. I|15l) 
becomes: 



4/(p) cos 



cosp 



2 
I 



A cosp + 2 cos(^) 



x Yj(p')co S p' = 6EWf(p). (16) 



Were it not for the A-term, the spectrum of eigenval- 
ues 5E( 2 >(Q,p) would be identical to that of two non- 
interacting soliton and anti-soliton with energy (2 cos q + 
2 cos q') = 4cos(Q/2) cosp with Q = q + q' and p = 
(q — q')/2. In the thermodynamic limit L — > oo at a 
given Q these form a continuum — 4cos(Q/2) < SE cont < 
4cos(Q/2). Hence, defining £ = SE& /4cos(Q/2), the 
bound states will be searched in the region \£\ > 1. Eq. 
lfl^|) can be solved formally as: 



M = r 



Ccosp + 1 



cosp — £ 
where C = A/[2cos(Q/2)] and 



1 

r = — J dp f(p) cosp 



(17) 



(18) 



Now, re- inserting l|17|) into l|18|) . self-consistency demands 
that either r = or: 

1 r , Ccosp +1 

1 = — / dp cosp — . 

27r cosp — t 

Evaluating the integral the energy equation becomes: 



1 C£ + l 
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\£\* 



C£ 



With A > one has a solution either for £ > 1 or for 
£ < min(— 1, — C _1 ). Under these conditions the solution 
reads £ = -(1 + C 2 )/2C. The condition £ < -C" 1 should 
be imposed for C < 1 but then the requirement would 
be C > 1. So we can only accept the case C > 1 for 
which £ < — 1 is always fullfilled. It can be also seen that 
Tlipfip) = is automatically satisfied if we build /(p) 
from eq. I|17|) ■ The restriction C > 1 can be re-expressed as 
cos(Q/2) < A/2 and defines a characteristic wavenumbcr 
Q* above which the bound state with energy 



SE^iQ) 



A -A C ° s2 lt 



> 



emerges from the continuum. This is precisely the result 
of eq. itTTl) . 
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